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Abstract 

We use ah initio total-energy calculations to predict the existence of polarons 
in semiconducting carbon nanotubes (CNTs). We find that the CNTs' band 
edge energies vary linearly and the elastic energy increases quadratically with 
both radial and with axial distortions, leading to the spontaneous formation of 
polarons. Using a continuum model parametrized by the ah initio calculations, 
we estimate electron and hole polaron lengths, energies and effective masses 
and analyze their complex dependence on CNT geometry. Implications of 
polaron effects on recently observed electro- and opto-mechanical behavior of 
CNTs are discussed. 

PACS numbers: 71.38. +i, 71.20.Tx, 71. 15. Mb 
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Typeset using REVT^ 



Carbon nanotubes (CNTs) have recently attracted a great deal of interest for their un- 
usual electronic and mechanical properties In this context, the interplay between me- 
chanical distortions and electronic structure plays a central role. Mechanical distortions 
can be generally classified as externally applied or spontaneous. Effects of externally ap- 
plied distortions such as twisting, bending and axial compression of CNTs on their electronic 
structure have been the subject of many studies Spontaneous distortions are usually 

related to strong electron-phonon interactions, and classic examples are polaron formation 
in ionic solids 0, Peierls distortion in ID metals [§| and related excitations (solitons and po- 
larons) in conjugated polymers 0. Symmetry-breaking distortions in fullerenes and CNTs 
have also been considered by several authors [|T0|, pJ] , p!^ , p!^ , p!^ ,p!5[] . 



In this work we show that an extra electron or hole in a CNT causes a completely 
different kind of spontaneous distortion: a combined radial (breathing-mode-like) and axial 
distortion. This perturbation causes the band edge energies to vary linearly and the elastic 
energy to increase quadratically with the distortion parameters. Therefore, the total energy 
of the system has a minimum at a nonzero value of the distortion parameter and a polaron 
is formed. 

The distortions considered here are changes in the CNT radius and length, characterized 
by the radial and axial strain parameters 

Sr = and = — , (1) 

respectively, where _Ro is the equilibrium radius for a neutral, undistorted CNT and io is 
the equilibrium length of its unit cell. Consider a semiconducting CNT with a single extra 
electron at the bottom of the conduction band. Allowing 6^,2 to depend on z, the CNT axis 



direction, we write the change in total energy caused by this extra electron as |T6| : 
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where V'(^) is the electronic wavefunction and frieff is the electron effective mass; kr^z and 
Xr^z are the effective spring constants per unit length and the electron-phonon coupling 
constants relative to the purely radial (r) and purely axial {z) strains; k^z is the spring 
constant relative to coupled radial and axial strains. Minimizing (§) with respect to ip* and 
z leads to the following expressions for the axial and radial strains: 



\z,rkrz Xr,zkz,r) 

and to the nonlinear Schrodinger equation 



ri^ = Cr,zri^ (3) 
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with (7 = (2meff/h'^){\rCr + XzCz)- Eq. (|) admits the following bound normalized solution: 

= sech(a2;), (5) 

where the inverse polaron length, a, and its binding energy, e eff, are given by 

C & 

The resulting maximum axial and radial distortions are e^^^ = aC^^zl'^- 

The polaron mass can be estimated using a semi-classical description of the electron 
motion. Assuming that the polaron propagates with velocity v^oi without changing its 
characteristic shape, we write the total kinetic energy of the electron-lattice system as the 
sum of the energy of an electron propagating in the conduction band with velocity v^oi plus 
the kinetic energy of the ions due to the propagation of the polaron: 

Both radial and axial components of the ionic velocities contribute to the ionic kinetic energy: 
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where Mj is the ionic mass and a is the hnear mass density of the CNT. Using Eqs. 
(H), we obtain the polaron mass: 



1 + 



2aaC^ 



15m, 



(8) 



and 



(9) 



Exactly the same results are obtained for an extra hole in the valence band, with a 
reversed sign in the definition of A, as discussed below. All quantities defined above are 
functions only of kr,z,rz, K,z, ^eff, Ro, and cr. We determine Ro, Eq and a from the 
CNT optimized geometry, and rrieff is taken from tight- binding (TB) calculations with 
a hopping parameter of 70 = —3.03 eV. 

We obtain kr,z,rz and Xr^z from ab initio calculations. Our calculations are performed 
using the SIESTA code, a numerical-atomic-orbital method based on density functional 
theory. This technique has been successfully applied to a number of studies involving nan- 
otubes P,p!B|. We use a generalized gradient approximation for exchange and correlation and 
norm conserving pseudopotentials. A split-valence double-C basis of pseudoatomic orbitals 
with an orbital confining energy of 0.3 eV and an energy cutoff of 300 Ry for the fast Fourier 
transform integration mesh are used. In order to probe the dependence of polaron properties 
on CNT geometry, we perform calculations for (11,0) and (7,0) zigzag nanotubes. These are 



within the range of recently reported experimental values for nanotube diameters The 
k-point sampling is composed of six k-points in the CNT axis direction, which allows us to 
use a minimal hexagonal supercell of 44 atoms for a (11,0) CNT and 28 atoms for a (7,0) 
CNT. The in-plane lattice parameter was chosen to be large enough (30 A) to ensure that 
there is negligible interaction between periodic CNT images. 

The calculations are performed for uniformly distorted neutral CNTs. Then e,. and 
become independent of z and the change in total energy per unit cell can be written as 
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AEtot = y-^K4 + + Kzere^j . (10) 

We vary the CNT radius while keeping = 0, and we vary the CNT unit cell length 
while keeping = 0: kr and kz are obtained, respectively, from the curvature on the AEtot 
vs. er,z plots, as shown in Fig. |l](a) for the (7,0) CNT and in Fig. | (b) for the (11,0) CNT. 
Then, by relaxing the unit cell length for a given nonzero value of the radial strain e^, kj-z is 
obtained. The resulting values for kr, kz and krz for both (11,0) and (7,0) CNTs are shown 
in Table |. 

The electron-phonon coupling constants for electrons (holes), A^'-^-', are obtained from 
the linear variation in the conduction (valence) band edge energies for a given strain energy 

€-r,z ■ 

AE,f ) = ±Kl^hr,z, (11) 

where the positive sign is for electrons and the negative sign is for holes. Plots for AE^^^^ 
for both (11,0) and (7,0) CNTs are shown in Figs, ^(a) and (b), respectively. Values for all 
electron-phonon coupling constants are displayed in Table |. 

The calculated polaron signatures (lengths, binding energies, masses and maximum dis- 
tortions) calculated from the ab initio elastic and electron-phonon constants are also pre- 
sented in Table |. The differences between the results for the (11,0) and (7,0) tubes allow 
us to anticipate a rich dependence on the CNT's diameter and chirality. This dependence 
can be understood as a superposition of two contributions: a "classical" and a "quantum" 
contribution. The "classical" contribution comes simply from the dependence of the elastic 
constants on the CNT's diameter. From Table |, one sees that all elastic constants are 
smaller for the thinner (7,0) CNT. This effect can be understood in simple terms by consid- 
ering CNTs with different diameters where the same percentual radial distortion is applied. 
Each bond in the CNTs' zigzag chains will be deformed by the same amount, but the thinner 
CNTs have less bonds along the chain. Considering the elastic energy to be proportional to 
the number of deformed bonds, the thinner CNTs will have smaller A;'s. A similar argument 
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applies to the axial elastic constants. So, considering this "classical" argument alone, one 
would expect stronger polaron signatures for thinner CNTs. 

The "quantum" contribution comes from the different possible signs and magnitudes of 
the electron-phonon coupling constants A^^^^. As one sees from Table I, it is difficult to 
identify simple trends of the A's with geometry. Simpler TB models predict a constant value 
of A for zig-zag tubes under uniaxial strain P|. This is probably an effect of curvature and 



rehybridization of cr — vr orbitals in small-diameter CNTs, which are generally neglected 
in TB calculations. 

Analysis of polaron signatures in Table | might suggest that the observation of polarons 
in CNTs would be very difficult. Although we have not attempted to reach particular 
combinations of chirality and carrier type to achieve the strongest polaron signatures, the 
studied examples yield binding energies at most ~ 0.3 meV, and the largest polaron mass 
~ 0.2% larger than the free electron mass for the hole polaron in the (7,0) CNT. 

However, these are signatures for a single polaron. The situation changes completely 
when we consider collective polaron effects. One example is the CNT length variation due 
to polarons. The total length variation on an CNT caused by the presence of a single polaron 



is: 



a6^^ = / dz e,{z) = C, (12) 
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For the hole polaron in the (7,0) CNT, Cz = 0.02 A. Therefore, a modest number of 500 
polarons would cause a sizeable 10 A variation in the nanotube length, large enough to be 
observed, for instance, in AFM or STM experiments where CNTs are used as probes. A 
related quantity is the strain- charge coefficient (SCC) pT| , 

Ai/i _NM^_NC, 

Ay - io " ^0 ' ^ ^ 

where Ai/i is the fractional change in length caused by a Ay change in the concentration 
of injected charge per carbon atom and is the number of atoms in the unit cell. Our 
calculated values for the SCC are also presented in Table |. They are comparable to the 
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experimentally measured values of 0.17 (for low charge injection) in the context of electro- 



mechanical actuation in CNTs This suggests that polaron formation may contribute 
significantly to the observed actuation, although it may not be the only driving mechanism, 
since experiments are usually undertaken in much more complex environments where inter- 
tube interactions and other collective effects may be important. 



Opto- mechanical effects in CNTs have also been observed |2^. Our results allow us to 
predict strong or weak axial distortions as an elastic response to light, depending on the 
signs of Cz for electrons and holes. Consider, for instance, electron-hole pairs generated by 
light in (11,0) and (7,0) CNTs. From the signs of Cz (the same signs of e^"^' in Table |), we 
can see that an electron polaron causes an axial expansion in both (11,0) and (7,0) CNTs. 
On the other hand, a hole polaron causes an expansion in the (11,0) and a contraction in 
the (7,0). Therefore, for the (11,0) CNT, the axial mechanical effects of the electron and 
hole will add up and the CNT will have a strong elastic response to light. On the other 
hand, electron and hole strains will partly cancel out in the (7,0) CNT, leading to weaker 
opto-mechanical effects. 

In conclusion, we predict the existence of polarons in semiconducting CNTs. Polaron 
properties are estimated from ah initio total-energy calculations for neutral uniformly dis- 
torted CNTs within a continuum approach. The complex dependence of polaron signatures 
on CNT chirality is understood as a combination of classical (elastic) and quantum effects. 
We show that collective polaron effects may have implications in the context of recently 
observed electro-mechanical |21| and opto-mechanical p2[ activities in CNTs. In particular, 
we predict the existence of two types of nanotubes regarding their elastic response (strong 
or weak) to light. A quantitative description of these effects should involve inclusion of 
collective behavior, inter-tube interactions and electrostatic interactions between electrons 
and holes should be taken into account into the model calculations. 
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TABLES 



TABLE I. Calculated quantities for a polaron in semiconducting zigzag CNTs. 



Parameter 


(11,0) CNT 


(7,0) CNT 




electron 


hole 


electron 


hole 


Xr (eV) 


0.99 


7.30 


-5.65 


-3.91 


A. (eV) 


-8.29 


-2.38 


-1.55 


7.76 


kr (eV/A) 


605 


605 


194 


194 


(eV/A) 


644 


644 


198 


198 


krz (CV/A) 


142 


142 


78 


78 


e (meV) 


-7.7 X 10-2 


-7.1 X 10-2 


-3.5x10-2 


-2.7 xlO-i 


L(A) 


390 


404 


586 


209 




1.00020 


1.00012 


1.00001 


1.00221 


(max) 


-6.291 X 10-6 


-1.684 X 10-5 


1.205 X 10-5 


3.413 X 10-5 


(max) 


1.783 X 10-^ 


8.212 X 10-6 


9.623 X 10-^ 


-5.275 X 10-5 


1 Aj/ 1 


0.142 


0.068 


0.074 


0.146 
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FIGURES 



(a) 



^ -4321.6 




(b) 
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Fig. 1 - Veriesimo Alves et al. 

FIG. 1. Variation in total energy as a function of (triangles) and (squares) for (a) the (7,0) 
and (b) the (11,0) CNTs. The solid curve is a 3^'^-order polynomial fit to the data. The values of 
kr,z,rz are obtained from the 2"*^ derivative of the polynomial at the minimum. 
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Fig 2 - Vens^imo Alve^ etnl 



FIG. 2. Band edge energies as a function of (squares) and (circles) for (a) the (7,0) and 
(b) the (11,0) CNT. Data around the regions er,z — are well fitted by straight lines. 
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